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Most d° transition metal (TM) oxides exhibit anomalously large Born dynamical charges associ- 
ated with off-centering or motion of atoms along the TM-O chains. To understand their chemical 
origin, we introduce "Wannier Orbital Overlap Population" (WOOP) and 'Wannier ond Orbital 
Position Population" (WOPP) based on the Wannier function based description of electronic struc- 
ture obtained within first-principles density functional theory. We apply these concepts in a precise 
analysis of anomalous dynamical charges in PbTiOs, BaTi03 and BaZrOs in the cubic perovskitc 
structure. Determining contributions of different atomic orbitals to the dynamical charge and their 
break-up into local polarizability, charge transfer and covalency, we find that p orbitals of oxygen 
perpendicular to the -TM-O- chain contribute most prominently to the anomalous charge, by facili- 
tating a transfer of tiny electronic charge through one unit cell from one TM atom to the next. Our 
results explain why the corner-shared linkage of TM06 octahedra, as in the perovskite structure, is 
ideal for large dynamical charges and hence for ferroelectricity. 

I. INTRODUCTION 

Transition metal oxides (TMOs) exhibit a rich variety of structural, electric and magnetic properties^, 
owing mostly to the d electronic states of the transition metal. Partially occupied d states of a TM 
in its oxide result in different possibilities of magnetic ordering and properties, where as the TMOs 
with unoccupied d states usually are non-magnetic insulators with structural and dielectric properties 
of technological importance 3 -. Most ferroelectric (ABO3) oxides in the perovskite structure^ have a d° 
transition metal at the i?-site, and interaction between d states of TM with p states of O are known 
to be crucial to their properties^. Many binary d° TM oxides such as Ti02 and Zr02 hold promise as 
high-fc dielectric materials^ essential to the future ultra-thin semiconductor devices. Understanding how 
precisely the d° state of a TM atom and its interaction with oxygen gives rise to these properties is 
fundamentally interesting and will also help in design of better materials. 

Dielectric properties are intimately related to electric polarization (dipole moment per unit volume) of 
insulators. It is spontaneously present in a ferroelectric material^ even in the absence of external fields, 
and its coupling with structure, electric and elastic fields are central to technologically important prop- 
erties of these materials. Dynamical charges, known as Born effective charges(BEC), yield polarization 
induced due to small displacements of atoms from their equilibrium positions. Anomalous charge is eval- 
uated as the difference between the dynamical charge associated with an atom and the nominal charge 
represented by it's oxidation state. They are also the coupling constants between atomic displacements 
and electric field, and are relevant to the properties of materials in two important ways: (a) they often 
are anomalously large and indicators^ of the ferroelectric structural instabilities or low frequency (soft) 
polar phonon modes and (b) static dielectric (piezoelectric) response depends quadratically (linearly) 
on the effective charges and is dominated by the softest phonons indicated by the anomalous effective 
charges^. 

First-principles calculations based on density functional theory (DFT) have been extensively used 
in determining dielectric^ and ferroelectric^, properties of the d° transition metal oxides. Within first- 
principles DFT framework, BECs are commonly calculated using two methods: (a) Berry phase based 
expression of polarization and (b) DFT-linear response. Both approaches are based in the A:— space and do 
not readily yield a picture in real-space of the flow of electronic charge in terms of bonding. Band-by-band 
decomposition 10 of effective charges using the linear response method was used to estimate contributions 
of different atomic orbitals to their anomalously large values, highlighting the role of hybridization between 
d states of TM with p states of oxygen. The band- by-band decomposition 11 of BECs associated with 
the Ba, Ti and O atoms in BaTi03, in groups of bands corresponding to different maximally localized 
Wannier functions^, was found to be in good agreement with that obtained using the linear response 
based approach by Ghosez etal^. Recently, analysis based on Wannier function o 12 ' 13 has indicated the 
important role played by the oxygen orbitals in giving rise to large anomalous Z* of the TM atoms in 
TMOs. However, a quantitative estimation of various contributions such as covalency, local polarizability 
and charge transfer to Z* is still lacking. 

In this work, we link Wannier function based description with atomic orbitals to quantify various 
chemical contributions to Z*. Polarization along x— direction is given by the trace of PxP, where x 



2 




FIG. 1: Two unit cells of a cubic ABO3 perovskite system. Labels in this figure are used in discussion. 

is the position operator and P an operator of projection into subspace of occupied states. A natural 
microscopic picture of its makeup can be obtained by diagonalization of the PxP operator. For systems 
periodic in 1-dimension, cigcnfunctions of PxP operator are the Wannier functions^ (or hermaphrodite 
orbitals of 3-D periodic systems^). For 3-dimensional periodic systems, PxP, PyP and PzP do not 
generally commute and can be diagonalized simultaneously only approximately, resulting in maximally 
localized Wannier functions^. In analogy with crystal orbital overlap population^, based on projec- 
tions of Wannier functions onto atomic orbitals, we propose here concepts of Wannier Orbital Overlap 
Population (WOOP) and Wannier Orbital Position Population (WOPP), which are used to determine 
quantitatively the contributions of different electronic mechanisms to anomalous dynamical charges. In 
principle, Mullikeni^ charge population becomes a subset of WOOP, if the Wannier functions are replaced 
by an atomic orbital basis set in which the electronic wave functions are expanded. On the other hand. 
WOPP introduces a new scheme to accurately estimate the degree of covalency. 

In section UH we describe the systems analyzed in this work. We review in section IIIII the definition 
of BEC, a method to construct localized Wannier Functions (WFs), and introduce WOOP and WOPP. 
In Section IIV( we present results of our analysis and break-up of the effective charge into different 
mechanisms, and finally conclude in section [Vj 

II. SYSTEMS 

PbTiOs (PTO) is an end member [x = 1) of PbZi'i-^Ti^Oa (PZT), the most commonly used techno- 
logical ferroelectric material. It occurs in the perovskite structure with a tetragonal ferroelectric phase 
up to a relatively high Curie temperature T c = 763 K, above which it transforms to the cubic paraelectric 
phase. While the stereochemical activity of lone pair of 6s electrons in Pb is known to be crucial for 
its properties, it shares with other perovskite ferroelectrics the d° TM-oxygen interaction contributing 
to ferroelectricity. It has been studied extensively using first-principles calculations^. As indicators of 
ferroelectric instability, effective charges of Pb and Ti in its cubic paraelectric phase (a — 3.97A) are 3.9 
and 7.1 respectively, which are analyzed here. 

BaTi03 (BTO), a related material, undergoes a sequence of structural transitions from cubic to tetrag- 
onal, then to orthorhombic and subsequently to rhombohcdral phase as temperature is reduced. In the 
cubic structure (a = 4.00A), BECs of Ba and Ti are 2.74 and 7.2 respectively. Polarization of BaTi03 is 
relatively small partly because Ba does not have any lone-pair of electrons and does not contribute much 
to ferroelectricity. A comparative study of PbTi03 and BaTi03 allows us to identify the special role of 
Pb at A-site. 

Finally, we analyze effective charges in BaZrOs (BZO), which remains stable in the paraelectric phase 
at all T, with the cubic perovskite structure (a = 4.19A). Effective charges of Ba and Zr are 2.74 and 5.7 
respectively. Comparison of various contributions to anomalous charges of Zr and Ti allows us to single 
out the chemical mechanism which leads to ferroelectricity in BaTiOs, but not in BaZr03. 

III. METHODS 

All calculations reported here have been carried out with ABINIT— implementation of DFT using 
Tctcr extended norm-conserving pseudopotentials and a plane wave basis. The s and p semi-core states 
of Ti, Zr and Ba are included in the valence. An energy cutoff of 80 Ry on the plane-wave basis is 
used in representing Kohn-Sham wave functions. Integrals over the Brillouin zone have been sampled 
with a 6x6x6 uniform mesh of k-points. Wannier functions have been constructed using post-processing 
subroutines interfaced with ABINIT code, developed by the authors^. 
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A. Born effective charge 



Elements of the BECs of an atom (I) are tensors, whose elements Z* a ^ 7 are the coupling constants 
between atomic displacements (phonon) and electric field. Force F felt by an atom / in the presence of 
electric field E is: 



Z a/3,I E 0, 



where as the dipole moment induced by displacement u of an atom is: 

Afi a = Z* apl u I(j . 



(1) 



(2) 



Corresponding change in the polarization is Afi/V, V being the volume of unit cell. In this work, we use 
the relation ^ in analyzing the BEC. 

A well localized WF W n (r) (constructed in the scheme described below) enables evaluation of expec- 
tation values of position operator: 



< x > 7 



xW*(r,0)W n (r,0)dr. 



(3) 



The total electronic contribution to the dipole moment per unit cell is obtained by summing over all the 
bands: fi x = ef ^2 n < x >„, e being the electronic charge and / the occupation number of the nth band. 
The BEC can be determined from the changes in the dipole moment resulting from the off-centering of 
an atom. For example, BEC associated with dipole induced along x direction by a small displacement 
u x of an ion / is given by: 



7* 



(4) 



where Zf n is the ionic charge (Zy n = number of valence electrons of the given atom treated with a 
pseudopotential) . A/i x is calculated from the WCs as: 



*l4 = e 



N 

zZ <x 



N 



< x > 



nd 



(5) 



where < x >f and < x >" d are the expectation values of x obtained using the i-th WF with and 
without the atomic displacement respectively. WFs thus allow decomposition of the dynamical charge 
into contributions from different bonding orbitals spanning the occupied subspace. 



B. Wannier functions 



Wannier functions(WF) are obtained as Fourier transform of Bloch functions of occupied bands. At 
any wave vector k, any unitary transformation of Bloch functions within the occupied band subspace 
is a valid description of electronic structure in the sense that the total energy is invariant under such a 
transformation. As a result, WFs are highly non-unique. Nevertheless, there are always some special 
choices of these unitary rotations, also known as gauge, which yield highly localized sets of WF o 16 ' 21 . 

In this work, we use the Wannier type localized orbitals (WLO) described in Re£i: 

Wj(r, R) = £ e 4k - (r ~ R) £ M 7 % k (r), (6) 
k i 

where the wave functions {u;k(r)}, obtained from a DFT calculation, are the cell periodic part of Bloch 
functions. Rotation of these functions by the unitary matrix M k (at each k) is aimed at obtaining a set 
of wave functions that are smooth and periodic in reciprocal space, a criterion necessary to obtain well 
localized WFs. 

In ionic systems, WLOs are typically centered on atoms, and their delocalized features extending 
to nearest neighbour atoms indicate the degree of covalency. In strongly covalent systems, WLOs are 
typically bond-centered, while the atom-centered WLOs can be obtained by expanding the subspace of 
electronic states to include the unoccupied, antibonding states. As demonstrated in Fig-[2] the Ol-centered 
[see Fig-[T] WLO with 2p y orbital character and the Ti-centered WLO with 3d xy orbital character, are 
constructed using an extended subspace above the Fermi level, in cubic PTO. The width or "fatness"— of 
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FIG. 2: Band structure of PbTi03 with (a) 2p y and (b) 3dxy characters centered at Ol and Ti respectively. 
The 2p y WLO at 01, constructed (c) within the occupied subspace, (d) in an extended subspace including bands 
above Ef. (e) The 3d xy WLO at Ti constructed in the extended subspace. 



the bands in Fig-^a) and (b) is linearly proportional to their contribution to the two WLOs respectively 
while Fig-[2jc) shows the 2p y orbital of 01 constructed within the occupied subspace. When constructed 
in the extended subspace, the antibonding contribution between the orbitals 2p y of 01 and the 3d xy of 
Ti from the bands above Fermi level in Fig-[2ja) destructively interfere with the corresponding bonding 
contributions from the bands below Fermi level to yield a pure 2p y like WLO shown in Fig-[5£d). Whereas, 
for the pure Ti 3d xy , shown in Fig-[2je), the majority contribution comes from the band above Fermi level, 
as clearly indicated in Fig-[^b). 

In the rest of this paper, we refer to the atom-centered WLOs constructed in the extended subspace 
as atomic orbitals (AO) and WLOs constructed within the occupied subspace as WFs. 



C. Wannier Orbital Overlap Population (WOOP) and Wannier Orbital Position Population 

(WOPP) 

A picture of bonding in a crystal obtained with Wannier functions is qualitative. To link it quantita- 
tively to different chemical mechanisms of bonding, and we now introduce tools based on their projections 
on atomic orbitals. For clarity, we denote indices of WFs and atomic orbitals (AO's) in upper and lower 
case letters respectively. 

In the first step, we express the WFs in terms of AOs {|</>j)} as: \Wi) = Ylj=i ^ij\4>j): where the AOs, 
N in number, may not necessarily form an orthonormal basis. Cjj's can be determined by solving N 
simultaneous linear equations of the form: 

N 

^Snjdj = (ct>m\Wl), 
i=i 

where S m j = (4> rn \4>j). Normalization of \Wi) in terms of AOs is: (Wi\Wi) = Y^hn^h^^Cim- 
Evidently, a choice of AOs is quite important in obtaining the normalization of WFs accurately. 

The part of the summand corresponding to I = m in the normalization of \Wi) gives the population of 
electrons at \cj>i) due to \ Wj). We generalize this to define "Wannier orbital overlap population" (WOOP) 
associated with the 7-th WF as: 

Bfm = CnSi rn Ci m . (8) 



(7) 
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Choosing I and to corresponding to all the AOs centered on an atom, say Al in Fig-[31 the population of 
electrons, which are centered on Al and contribute to the Wi, can be evaluated as: 

QAl= E B L- (9) 
ImeAl 

Subsequently, the number of electrons populated in atom Al can be estimated as Y] T Qai- 

Center of \Wi) can now be expressed in terms of the AOs: (x)i = J2 lm C* n Xi m Ci m , where Xi m — 
{4>i\ x \<t>m) ■ We refer to the summand as "Wannier orbital position population" (WOPP) and denote it as: 

D lm = C* n Xi rn C Im . (10) 

This leads to the net contribution to Z* from \Wi) as: 

N 

Z *=7T E (Dt-D?* 1 )- (11) 

l,m—l 

Through appropriate choices of I and to in ecin- lfTTj) . the Z* 1 can be easily decomposed into contributions 
based on different chemical mechanisms responsible. These are classified into five broad categories: rigid 
shift (RS), charge transfer (CT)i£, local polarizability (LP), covalency (COV) and other nonlocal changes 
(NC). For the WF centered on Al (see schematic in Fig-[3]), these are defined by the following conditions 
on I and to: 

1. RS: I = to, {lm} G Ai, AQ\ { = 0, i = 1, 2, 3; 

2. CT: I = m, {lm} e Ai, AQ^ ^ 0, i = 1, 2, 3; 

3. LP: I ^ to, {/to} e A; 

4. COV: / ^ to, I e A, to ^ A or to e A, / ^ A; 

5. NC: I ^ m,{lm} £ A. (12) 

RS and CT are contributed from the diagonal blocks in Fig. [3] (identical AOs) and can be differentiated 
based on the change in the population (AQ), obtained through WOOP. Whereas, Z* from the LP arises 
solely from the coupling between the AOs with different parities centered on the same atom, for which 
the inter orbital transition dipole moment is nonzero. Due to orthonormality of AOs centered on the 
same atom, contribution of LP to Z* is not effected by the rigid translation of AOs. 

Contributions to Z* from COV and NC are both due to transition dipole moments between orbitals 
centered on different atoms. In Fig-[3] where |uT) has highest projection with AOs of Al, COV gives 
covalent contribution to Z* due to transition dipole moments between the AOs of Al and it's nearest 
neighbours A2 and A3. COV contributions can be both bonding as well as antibonding type in nature. 
Finite COV as well as CT contributions to BEC, involving same set of AOs of two neighbouring atoms, 
can be indicative of a possible conjugation scenario, as shown in section- ITVl NC gives the contribution 
to Z* from covalency associated with transition dipole moments between the AOs of A2 and A3[Fig-[3]. 
In general, NC contribution to Z* is much smaller than the other four categories as none of the AOs 
involved in NC belong to Al. 




FIG. 3: A schematic grouping of / and m to decompose Z* due to Wi(RS: Rigid shift, CT: Charge transfer, LP: 
Local polarizability, COV: Covalency, NC: Nonlocal changes). 
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IV. RESULTS AND DISCUSSION 

We first present WOPP based analysis of Zt ri in PTO and compare it with that of Zt ri in BTO. 
Next, we analyze Z* ph of PTO and bring out the special role played by Pb through comparison with Ba. 
Finally, we present a comparative analysis of BTO and BZO to identify factors crucial to ferroelectricity 
in BTO and to the lack of it in BZO. In all these cases, cations are off-centered (with respect to the cubic 
structure) along x— axis by 1 % of the lattice constant. 

In addition to a nominal charge of 4 a.u., the anomalous part of Z1^ i in PbTi03 is about 3.1 a.u. To 
this anomalous part, a positive contribution of 4.2 a.u. is from the orbitals of 01, which belongs to the 
Ti-0 chain along the x-axis, and a negative contribution of -1.1 a.u. from the orbitals of 02 and 03 
in the Ti-0 plane perpendicular to the x-axis. Almost 70% of the total positive contribution (4.2 a.u.) 
is made by the 2p y and 2p z orbitals of 01 (both are equal by symmetry), which are perpendicular to 
the Ti-0 chain along x-axis. The Ol-centered WFs with 2p y orbital character are shown in Fig-HJ The 
majority of their contributions (see Table-fl} to Z1 Ti is due to inter-atomic charge transfer across the unit 
cell from the 3d xy and 3d xz orbitals of Til to those of Ti2. WOOP analysis for these WFs clearly shows 
a decrease and increase in electron population at Til and Ti2 respectively, with negligible change in the 
population of orbitals at 01. 

Bonding in PTO is certainly not purely ionic: a small but non-zero covalent contribution of about 
0.3 a.u. to Z* from mixing between the 01 and Ti orbitals indicates a weak tt like bonding along the 
Til-01-Ti2 chain which, is strengthened in the ferroelectric phase. The rest 30% of the total positive 
contribution to the anomalous part of Zt ri is due to the 2s ± 2px hybridized orbitals of 01. Unlike the 
2p y and 2p z orbitals which are perpendicular to the -O-Ti-O- chain direction, covalent contributions from 
these orbitals arises from a type bonding between Ti and O and are comparable to that from charge 
transfer. These WFs (shown in Fig-[S]) are centered between O and Ti atoms: the 2s + 2p x orbital is 
peaked in the bond path between 01 and Til where as 2s — 2p x peaked is peaked in the bond path 
between 01 and Ti2. The former contracts and the latter gets stretched in the ferroelectric phase. 
Consequently, in the ferroelectric phase the covalency between 01 and Til gets stronger at the cost of 
the same between 01 and Ti2. This is similar to ir conjugation found in long chain conducting polymers. 
As shown in Table-|ll this results into a larger covalent type contribution to Z^ i from 2s + 2p x than 
that from 2s — 2p x WF. These orbitals too facilitate the transfer of small electronic charge from Ti2 to 
Til, resulting in charge transfer contributions to Z^ i that are about the same as their contribution with 
covalency. Small contributions from local polarizability of oxygen to Z1 Ti from these orbitals correlate 
with the intra-atomic change in hybridization among the s,p and d orbitals of oxygen without altering 
the total charge population or WOOP at 01. Similarly, local polarizability of 02 (and changes in its 
hybridization) yield a small negative contributions to Z^. 

Interestingly, charge transfer between TM atoms at -B-site through O seems to be indirectly influenced 
by the atom A (or its size). Comparison of contributions to Zt ri in PTO and BTO (shown in Tables-fl] 
and HI) indicates that charge transfer from Ti2 to Til is higher in BaTi03 than in PbTi03. In contrast, 
the contributions from the 2s±2p x hybridized orbitals (a bonding) of 01 in BaTiOs are about the same 
as in PbTi03. As we will see below, the former difference is due to the fact that Pb also competes in 
contributing to the charge transfer through the same 01 orbitals. 



TABLE I: Various chemical contributions to and Z Pb in PbTi03.(CP:Charge transfer; RS:Rigid shift; LPLocal 
polarizability; COV:Covalency; NC:Non-local changes) 



WF 


CT 


LP 


COV NC 


Tot 






Decomposition of Z\ i 


from orbitals of 01 




2s + 2p x 


0.21 


-0.10 


0.41 0.01 


0.53 


2s - 2p x 


0.33 


-0.06 


0.30 0.01 


0.58 


2p y ,2p z 


1.10 


0.00 


0.27 0.06 


1.43 






Decomposition of Z Pb 


from orbitals of 02 




2p x 


0.12 


-0.08 


0.28 0.02 


0.34 


2p y ± 2s 


0.02 


-0.04 


0.03 -0.00 


0.01 


2p z 


0.38 


-0.03 


0.07 0.01 


0.43 






Decomposition of Z Pb 


from orbitals of Pb 




WF 


RS 


LP 


COV NC 


Tot 


6s 


-1.97 


0.58 


-0.02 -0.02 


-1.43 
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Anomalous part of Z* pb in PbTiC>3 is about 1.9 a.u., over a nominal charge of 2.0 a.u. Almost 25% 
of this anomalous part is a result of local polarizability at Pb due to emergence of moderate Pb 6p 
orbital character in occupied states of the ferroelectric phase, which is accompanied by a small decrease 
in occupancy of the 6s orbital. This is evident from the nature of deformation of the Pb-centered WF of 
s symmetry [see Fig-[6ta)] in the ferroelectric phase. The rest of the anomalous part is contributed by the 
two oxygens 02 and 03 in the (100) Ti-0 plane. 2p x and 2p z symmetry WFs centered on both of these 
oxygen atoms give positive contributions to Z Pb . Interaction between Pb and 02 orbitals is evident in 
the long-range features of the 2p x and 2p z like WF centered on 02 (shown in Fig-[^c) and (b)) in the 
ferroelectric phase. The contribution from 2p x is primarily covalent in nature involving a combination of 
6p orbitals of Pb, while the same from 2p z is mainly due to charge transfer from similar orbitals of Pb2 
to that of Pbl (see Table-|T|. 

In both BaTi03 and BaZrOs, the anomalous part of Z Ba is about 0.75 a.u. The hd orbitals of Ba, 
which are unoccupied in the paraelectric phase, mix with occupied states in the ferroelectric phase giving 
a positive contribution to Z Ba . These orbitals are relatively more localized than the 6s or 6p orbitals of 
Pb and hence less favourable for inter-atomic orbital interactions. Consequently, reduction in both charge 
transfer as well as covalent contributions involving the 2p orbitals of 02 and 03 and the hd orbitals of 
Ba is expected. Similarly, the local polarizability contribution from Ba due to hybridization between its 
hp and 5c? orbitals in the ferroelectric phase is also weaker. These are evident from comparison between 
contributions to Z* Ba (Tables-HQ and HTT|) and those to Z* pb (Tablc-|J). Our analysis of Z Zr (Table-HED) 
reveals two interesting trends: (a) contribution from covalency between Zr and O (c-type bonding) to 
anomalous part of Z* Zr is greater than that from covalency between Ti and O to and (b) charge 
transfer is much more significant in giving anomalous Z^, i than in the case of Z Zr . 

We sum up contributions from all the WFs and summarize the decomposition of anomalous Z*'s of 
A and B cations into various chemical mechanisms (see Tablc- ITV| . Comparison between BZO and BTO 
clearly reveals the role of long-range charge transfer between TM cations via p orbitals of oxygen (see 
Fig-[7]). This is analogous to super-exchange in magnetic oxide materials. As expected from the fact that 
the radius of Zr 4+ is larger than that of Ti 4+ , contribution to Z* from covalency is greater for Zr than for 
Ti. In case of Z Pb , the net charge transfer through the orbitals of 0(2,3) is the prime contributor, followed 
by the net contribution from covalency between the orbitals of 0(2,3) and Pb. The local polarizability of 
the 6s orbital of Pb also contributes sizeably to Z Pb , but is partly compensated by the local polarizability 
of oxygen. The cause for difference between Z Pb and Z Ba is most evident in charge transfer contribution, 
and to a lesser extent in contributions from their local polarizability and covalency with oxygen. 

V. CONCLUSION 

We presented a simple scheme to quantify different chemical bonding mechanisms through projections 
of Wannier functions on to atomic orbitals. This is particularly useful in separating the total BEC into 
various chemical contributions such as local polarizability, covalency and charge transfer. We find that 
anomalously large effective charge of Ti in BaTi03 and PbTi03 arises primarily from a dynamical transfer 




FIG. 4: Ol-centered WFs with symmetry of its 2p y orbital in cubic (a) BaTiOs and (b) BaZrOs (l)with and (2) 
without Ti/Zr displacement. The isosurfaces shown here, are for |W(r)| = 0.0045, which is typically about 10% 
of its maximum. 
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of tiny electronic charge through one unit cell from d orbitals of one Ti atom to its nearest neighbour Ti 
atom in the direction opposite to Ti displacement and secondarily from covalency (the latter is less than 
half of the former). In contrast, the anomalous part of effective charge of Pb has comparable contributions 
from charge transfer and covalency, of which the former is lacking in the case of Ba. 

We have identified the dynamical charge transfer between cations (see Fig. -[7]) (TM as well as the one 
at A-site) as the key factor for their potential to drive ferroelectricity. As this charge transfer occurs 
along -TM-O- chains through the p orbitals of oxygen that are perpendicular to the chains, it is clear that 
perovskite is a structure very special for ferroelectricity for two reasons: (a) it has -TM-O- chains running 
in all the three directions and (b) its TMOg octahedra are all corner-shared with two-fold coordination of 
oxygen atoms giving two of the p orbitals of oxygen for facilitating this charge transfer. More generally, 
the tools of WOOP and WOPP developed here, will be useful in characterization of chemical bonding in 
variety of difference systems and phenomena. 
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displacement in cubic PbTiOa. The isosurfaces are for |W(r)| = 0.0045, which is about 10% of its maximum. 
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FIG. 6: WFs with symmetry of (a) 6s orbital of Pb, (b) 2p z and (c) 2p x orbitals of Ol in cubic PbTiOa. Isosurfaces 
shown in (b) and (c) are for |W(r)| = 0.001. 
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TABLE II: Various chemical contributions to Z^i and Z* Ba in BaTiOa. 



WF 


CT 


LP 


COV NC 


Tot 






Decomposition of T^ ri 


from orbitals of Ol 




2s + 2p x 


0.19 


-0.08 


0.42 -0.01 


0.52 


2s - 2p x 


0.34 


-0.05 


0.29 -0.00 


0.58 


2p y ,2p z 


1.18 


-0.01 


0.26 0.06 


1.49 






Decomposition of Z* Ba 


from orbitals of 02 




2p x 


-0.01 


-0.05 


0.31 0.01 


0.26 


2s ± 2p y 


0.04 


-0.04 


0.06 -0.00 


0.06 


2p z 


0.11 


0.00 


0.11 -0.00 


0.22 






Decomposition of Z* Ba 


from orbitals of Ba 




WF 


RS 


LP 


COV NC 


Tot 


5p x 


-1.94 


0.08 


-0.02 0.00 


-1.88 


5p y ,5p z 


-1.94 


0.09 


-0.00 0.00 


-1.85 



FIG. 7: A schematic picture of interatomic charge transfer across one unit cell. 



TABLE IV: Comparison of various contributions to anomalous effective charges. 





Charge 


Covalency 


Local 






transfer 




polarizability 








(from orbitals of Ol) 




Ti (BTO) 


2.89 


1.23 


-0.15 




Zr (BZO) 


1.61 


1.57 


-0.28 








(02+03+A) 


(02+03) 


(A) 


Ba (BTO) 


0.36 


1.06 


-0.26 


0.26 


Pb (PTO) 


1.08 


0.80 


-0.38 


0.58 



TABLE III: Various chemical contributions to Z* Zr and Z* Ba in BaZr03. 



WF 


CT 


LP 


COV NC 


Tot 






Decomposition of Z* Zr 


from orbitals of Ol 




2s + 2p x 


-0.16 


-0.10 


0.63 0.00 


0.37 


2s - 2p x 


-0.01 


-0.10 


0.50 0.01 


0.40 


2p y ,2p z 


0.89 


-0.04 


0.22 0.03 


1.10 






Decomposition of Z* Ba 


from orbitals of 02 




2p x 


0.03 


-0.03 


0.23 0.02 


0.25 


2s ± 2p y 


0.04 


-0.03 


0.04 -0.00 


0.05 


2p z 


0.13 


0.00 


0.07 -0.01 


0.19 






Decomposition of Z* Ba 


from orbitals of Ba 




WF 


RS 


LP 


COV NC 


Tot 


5p x 


-1.93 


0.07 


-0.03 -0.00 


-1.89 


5p y ,5p z 


-1.94 


0.07 


-0.00 -0.00 


-1.87 



